The superfluidity mechanism of He II 
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Abstract: Based on a first principles treatment of the excitation states we 
show that superfluidity of superfluid 4 He (He II) results from a reduction in 
^vq , the number of phonon wavevector K states A/2 (K) to a level that is negligibly 

low when the fluid is confined e.g. in a narrow channel, yet wider than the 
helium atom correlation length, A. This is as a result of the K discretization, 
a manifestation of the quantum confinement effect (QCE). The predicted rela- 
tive viscosity of a confined superfluid has the characteristic order of magnitude 
of experimental data (< 10~ 6 ). Furthermore, we show that at the edges of 
the resulting energy gaps, the A/2(A') presents discontinuity. When its corre- 
sponding energy exceeds the (first) gap, the superfluid flow exhibits a critical 
velocity v c . Our evaluation of v c (d) versus the channel width d, constrained to 
■ satisfy energy conservation, is in good quantitative agreement with experimen- 

ted ' tal data for channels with d > 10 -6 m. Meanwhile, a sharp turn about K oc 

ff) , v c in A/2 (K) resembles very well that of the experimental overshoot data. For 

• narrower channels of d < 10 -6 m < A in which the phonon excitation picture 

<^^ ) becomes inadequate, we instead represent the excitation in terms of single 

atoms with an effective mass, which yields a v c (d) in close agreement with 
experiment. Accordingly, the reduction in the number of atomic states results 
in superfluidity. The theoretical finding in this work, which can be termed the 

S QCE superfluidity mechanism, provides a consistent explanation for this 
puzzling phenomenon, the non-dissipative, superfluidity motion, of He II and 
could have a significant impact also on the understanding of other superfluids. 
PACS numbers: 67, 67.20.+k, 67.90.+Z, 67.40, 67.40. Db, 67.40. Hf. 



I. INTRODUCTION 

It is one of the most puzzling natural phenomena that in the superfluid phase, He II, and 
in a narrow channel and below the critical velocity, liquid 4 He flows as if there is "no 
friction". Systematic experimental investigations have shown[ that the viscosity of 
He II is 10 6 ~ 10 11 times smaller than that of normal fluids. Although the phenomenon 
of superfluidity has been the subject of many theoretical studies for the past 70 years, 
including the well know work by London[||], Landau [ || and Feynman[ [jj, there is still not 
a satisfactory solution to the problem. In this work, we derive a new theoretical model for 
the superfluidity of He II based on the relevant experimental indications. 
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II. THERMAL EXCITATIONS IN THE SUPERFLUID 



To understand the superfluidity, one first of all can learn from the thermal excitation prop- 
erty of the superfluid, for dissipation is caused by thermal excitations. Experiments, in 
particular inelastic neutron scattering measurements [ ^-[^j, have provided basic evidence 
concerning the nature of the thermal excitations in superfluid He II. Three aspects, in partic- 
ular, are relevant here. First, the neutron scattering function S(q, u) of a (pure) superfluid 
of size ranging from the macroscopic [ to the nano scale [ is, at a given momen- 

tum transfer q up to at least 1.93 A , sharply peaked at hu>(q), Tiuj{q) being the neutron 
energy transfer. For q < 0.6 A -1 , the excitations in a macroscopic superfluid are known[ [l6| 

to be well described as elastic sound waves; here the "macroscopic" or "bulk", which will be 
assumed for the superfluid till Sec. refers to a dimension larger compared to the helium 
atom correlation length. The combined information above indicates that the sound waves 
in the He II superfluid are quantized; namely in terms of phonons, each having an energy 
quantum £(K) = Hlu(K), where Ti is the barred Planck constant, and u>(K) and K (= q) are 
the phonon frequency and scalar wavevector for the isotropic superfluid. Second, u)(K) has 
been accurately determined from neutron scattering experiments [ 0|; in the low K ( < 0.6 

A -1 ) region of concern here: 

S(K) = hcj(K) = hciK, (1) 
ci being the velocity of first sound of He II. Furthermore, thermal conductivity measurement [ 



17 has shown that in the superfluid at ~ 1 °K and in the small geometries in question, 
phonon-phonon collisions are basically absent and the phonons principally collide with the 
specimen boundaries. We can thus effectively represent the propagating phonon wave here 
by a free "particle" of mass m p h in real space, for which the de Broglie relation states 
£ = (TiK) 2 /m p h- Combining this with Eq. ([!]), we get m p h = TiKjc\. The Schrodinger 

equation for the phonon particle, e.g. in the x direction, is — = £(K x )u(x); its 

energy is obtained from Eq. (|l|). The amplitude solution consists of the planewaves 

u x {K x )=u e iK ". (2) 

Equations ([[])- (|J), obtained above based directly on experimental observations, give a first 
principles representation of the excitation states in a macroscopic superfluid sample, and 
will be the starting point for our subsequent investigation of the confined superfluid. For 
information only, we mention that the above solutions have been obtained in our microscopic 
theory of the superfluid of He II [ |l8| through solving the equation of motion of the superfluid 
atoms (the helium atomic wavepackets). Third, the number of phonon excitations (which 
can only have occurred between the discrete energy levels of the superfluid) increases in 
proportion to the concentration of the superfluid in He II and becomes 100% at zero °K, 
implying that the superfluid in thermal equilibrium is prominently excited, via phonons. 
This suggests that insofar as the superfluid motion is concerned, the phonon excitation 
cannot be arbitrarily omitted and, in arriving at a mechanism of superfluidity, we are obliged 
to take its effect into full account. 



III. THE THEORETICAL BASIS OF SUPERFLUIDITY MOTION 



We consider a superfluid under the condition of thermal equilibrium in steady flow motion 
along the z axis in a channel of width d and length L. Unless specified otherwise, we assume 
d « L; this represents a two-dimensional spatial confinement, 2D SC. We also address the 
cases of one-dimensional and three-dimensional spatial confinements, ID SC and 3D SC. For 
simplicity, we assume that T is low ( < 1 °K) so that He II is virtually a pure superfluid. 

In this and the next section we discuss wider channels and in Sec. [yjthe narrower channels, 
as the superfluidity mechanism appears to differ in detail for these two cases. 
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A. Phonon states in a confined superfluid below the first excitation gap 



In a superfluid confined to a channel (2D SC) with a d of macroscopic scale, by Sec. || the 
primary thermal excitation consists in phonons. It is readily seen that for this confined fluid 
(simplified to have a squared cross section) the phonon waves, Eq. (|^), are subject to the 
non-trivial lateral boundary conditions e lKxX \ x -Q — e lK * x \ x= d and similarly for y, and along 
the z axis the trivial one e lKzZ \ z= o = e %K * z \ z= Li as a result of constructive interference. The 
phonon wavevector ~K(K X , K yi K z ) in the x and z directions are therefore restricted to the 
discrete values JC n (K X i,K y j,K z i), with 

K xi = iK x , K yj = jK y , i,j = 0, ±1, . . . , ± (jV/3 ~ 1) ( 3 ) 



and K zl = lK z , I = 0, ±1, . . . , ± (JV * ^ , Kb) 



respectively; K n = |K„| = . + K* + K*,. Where, 



yj zV 

K x = K y = — , and (4) 



are the smallest K values permitted, the K quanta, in respective direction. Np = —, 
= x, y and N z = -, a being the He II interatomic spacing. The corresponding minimum 
excitation energies, or energy gaps, are A p h/h — uj x = u)(K x t i+ u) — ui(K X i) and uj z = 
w(K z n + x}) — uj(K z i). For the 2D SC, the K xi and K y j states are degenerate. For a ID SC, 
the K y j and the K z i states will be degenerate. For the low energy excitations of Eq. ([!]), 
we have 

A ph = hciK x = and (5) 



hc\ 

17' 



Suppose that d is small such that given an external perturbation of energy e s , 



£., < e. 



pi 



, = E A ph- ( 6 ) 



where £p h refers to the minimal excitation energy per atom permitted by the conditions (||) ; 
J2 a wm be specified in Sec. 1IIB. Then, no excitations can be caused along the K x and K y 



axes. On the other hand, since L is large, ui z is therefore effectively zero and the low K z i 
modes can undergo continuous excitations. 

For the oj(K) of Eq. ([[]), the number of states in an interval [u>(K n ), u>(K n + 8K)) is 
obtained by enumerating the number of J£.' n (K' xi , K' y j, K' zl ) points contained in the spherical 
shell K n < K' n < (K n + SK). Now for the minimum sphere with respect to the K X i modes 
which has the radius 

K 1 = K xX = 1-K x = K z l, 
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FIG. 1. Configuration of the phonon wavevector K states, shown for K < 2K\ = 2K X in the 
first quadrant in the K y — plane, for a superfluid conditioned as (a) or (b). (a). The superfiuid 
is 2D confined in a channel of dimensions d' 2 x L, with l/(d/L) = 8. Within K < Kx, K x being 
the minimal K x i value, the only permitted K states are those along the K z axis indicated by the 
triangles. After K > K x , an additional number of states are permitted (the triangles) along the 
K X i = K x grid line. The additional states in the first shell K x < K < K x + 1 • K z , indicated by 
the triangles on the AB line in the upper-right graph, give rise to a discontinuous increment in the 
A/2 (K) (cf. the inset of Fig. ^b). (b). The superfluid is in a large container of the dimensions L 3 , 
being thus un-confined. The permitted K states are indicated by the filled circles in the regions 
which are unoccupied by the K states of the 2D SC sample of (a), and by the triangles along the 
grid lines permitted for those of (a). If l/(d/L) — 5 x 10 4 as is usually employed in experiments, 
the number of additional K states - the filled circles - of the unconfined sample would be about 
10 7 times more densely populated than illustrated in the figure. 

the only permitted states in it are the K ra (0, 0, K z \) points on the K z axis (triangles below 
K x in Fig. |l]a). Then, accumulating these points in the above manner up to a given K 
value: K n — K z i — IK Z (< K\), here n = /, gives the total number of phonon wavevector 
states, or integrated density of states for the 2D SC: 

f £{'=-i*XM' = 3£ + 1=21 + 1, K<K X 

W0= E^ /d W+2E,=_ M ^'oo (7) 
I = (f + 1) + 4i| i= i, K = K 1 . 

Where, Sqqi stands for Siji\i=j = o; Siji equals 1 at each crossing point of the grid lines K X i, 
Kyj and K z i, and zero elsewhere. The density of states V 2 {uj) (= dA ^"^"^ = 9 a ^ n ^ ) 
is thus 

V 3 (u)/(dn/du) = { 2 £ K = Kt'. ^ 
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where dn/duj = L/2ttci = 1/K Z C\. The A/2 (if) states, which can be continuously excited, 
will strictly speaking dissipate the flow motion in question, however, only to an extremely 
small extent. One way of seeing the small extent is by comparing with a large (unconfined) 
sample of the dimensions L 3 . For the latter, the allowed phonon states are similarly defined 
by boundary conditions as Eqs. (^) except that d is replaced by L. Accordingly the if 
states for it are defined by the crossing points (filled circles in Fig. |l]b) of the grid lines 
(K x i, K y j, K z i), equally spaced by a K z in all axial directions, the volume quantum being 
if J. Within a sphere of radius if = nK z , of a volume f2 = ^f-K 3 , the total number of states 
and density of states for the unconfined sample are 

WO = ;§ = !»», withAVi^f (^) 3 ; (9) 

and V (w)/(dn/du) = inn 2 . (10) 

Comparing A/^ifi) of Eq. (Q) with the A/o(ifi) above, we obtain that up to if = K\ 
the number of states for the 2D confined sample is reduced by a relative amount (noting 
L/d»l, 2): 

A/~ 2 (ifi) „ 3 f 

As can be seen, r2,o (solid line in Fig. ||a) is very small for L/d >> f . For Lid ~ 5 x f0 4 , 
which is a typical sample dimension ratio used in superfluidity experiments [0-0], the number 
of phonon states is reduced by a factor r2.o ~ 10 -10 , to 10 10 times smaller. As we also see, the 
extremely small A/2 (if) corresponds to the removal of those states (filled circles in Fig. |l|b) 
which are permitted for the unconfined sample, lying in the interim spanned by the discrete 
states K X i, K y j of the 2D SC sample. It is noteworthy that the degeneracy of each discretized 
state is not increased. The extraordinary reduction in A2(if) in the reciprocal if space is 
also as expected directly from its relationship with the real space. We can readily inspect 
this relation by multiplying the maximum numbers of modes in all three axial directions 
given in Eqs. (||); we get the total number of states: Mg{K max ) = N x N y N z = ^f- = N s e, 
N s $ being the number of atoms in a volume Vg under 6D SC. This equality says that, 
reducing Vq = I? to V2 = d 2 L corresponds to removing the atom population by the factor 
jf^ = p| = (-j) = - the same conclusion as Eq. (|TT|) : the simple evaluation above 
shows in addition that A/2 (if) is reduced over the entire if volume. 



B. The superfluidity mechanism 

Consider a superfluid bulk of N s atoms and total mass M s = mN s , m being the 4 He 
atomic mass, set to flow motion at velocity v s . Its translation momentum and energy are 
thus P s = M s v s and E s = ^M s v 2 . Suppose now there presents a flow-wall interaction, 
or a viscous force. The flow energy will then be gradually dissipated and converted to an 
additional, disordered thermal vibration energy of the atoms, which by Sec. || are in the 
form of phonons, in the interfacial layers. The total flow-wall interaction consists of the 
collision of all of the interfacial atoms from the superfluid side, denoting their total number 
by Na, with all of those on the wall side. The flow energy per collision atom, as the e s in 
Eq. (||), is thus, 

E 1 

£s = it = 2^ (12) 

This physically represents that the superfluid atoms, being strongly correlated, will collide 
with the wall as if each has an effective mass: 



N A 



m. (13) 
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FIG. 2. (a). Ratio of the integrated phonon destiny of states of a (9D confined superfluid, 
A/9 (If 1), over that of an unconfined superfluid A/"o(ifi), as plotted as a function of L/d; d and L are 
the smallest and largest dimensions of the confined sample and L 3 of those of the unconfined. 0—1 
or 2 refers to a ID or 2D SC. (b). The theoretical result for A/a (if) v.s. / plotted up to a K — 2K X 
and hence I = 2K X /K Z = 2L/d, with d/L = 5 x 10 4 ; / = L/d and 2L/d are the first and second 
lowest K xi states. A/a(-?0 has discontinuous increments at iK x ,i = 1, 2, . . ., the discontinuity at K 1 
is shown on an enlarged scale in the inset. The circles indicate the experimental overshoot data 
(from Ref p^), O cxp , versus flow velocity v s . c ^ p . O CX p is scaled to O p i ot = 9 - 9 * lc|4 Q cxp j n the graph, 
so as to fall within the magnitude of A/2 (If). Meanwhile, v 3 , cxp is scaled to 5x 1 ^° v s .cx P = / so that 
the turning point of the experimental curve (i.e. v s ~ 15 = cm/sec «£ xp ) occurs at I = L/d — 5 x 10 4 . 
Despite of the arbitrary scaling, the valid information is that the variation of the theoretical A2(Z) 
v.s. / at about / = L/d closely resembles that of the experimental data: both C x P and AT-2(K) 
exhibit a smooth increase before / reaching the critical velocity (corresponding to I — L/d) and 
after that a much more rapid increase, (c). The corresponding phonon density of states T>e(u), 
divided by a constant N s L/2nci, for the 2D confined superfluid (6 — 2, solid curves and circles) 
and an unconfined superfluid (0 = 0, dashed curve). The dotted thin lines are guide for the eyes. 
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On the other hand, since the inert character of the helium atoms upon entering the He II 
state should not change with respect to the wall, we can regard that the flow atoms (with 
an effective mass), on passing a wall atom, do not interact with it except at the moment of 
collision. Then, the maximum excitation energy, tu(K(v s )), will be directly proportional to 



Vg for the flow-wall inelastic collision here; see further Sec. [IV B| . A successful excitation 
is dictated by the following factors. 1). The number of available states for excitation at 
each accessible K value, which is represented by T>g(w{K)). When its e s satisfies (||), the 
flow momentum, being along the z direction, will directly convert to the only allowed K z 
phonon momentum. On the other hand, when e s is able to excite the K z i = iK x modes 
(i 7^ Q), the excess phonons, produced first in the z direction, will cause a momentary stress 
difference between the z and the lateral (x, y) directions, and hence a momentum flow from 
the former to the latter. This will continue until, presumably within a relaxation period, 
the excess momentum is uniformly distributed amongst all the degenerate states (within a 
quantum uncertainty range), as dictated by the isotropy of the fluid. 2). The probability 
of excitation of a phonon of to(K) on perturbation. This is determined by the Boltzmann 
factor f(u),T) = exp(— %uj /ksT), assuming a one-phonon process for the low temperature 
involved. The sum of f{u,T) over all excitation events in a certain time interval At should 
equal [< n > (T + AT)- < n > (T)], < n > (T) being the Planck distribution at T, and 
AT the temperature increment due to the excess phonons in At. 3). The rate of head- 
on collisions in which a flow atom loses its entire translation energy. We define 7 as the 
number of head-on collisions an interfacial atom commits per unit travelling-distance, and 
hence T = j£ that after travelling a distance £. 4). There present in the interfacial layer 
Na simultaneous collisions at any one time. Multiplying the above factors and summing 
over all accessible states, we then obtain the total excess excitation energy produced by the 
flow-wall interaction over a travelling distance £ in a channel of #D SC: 

E ph .e(Lu(K))=TN A . e £ hu'T>g(uj')f(ij',T) 

lu'=lu(K')<w(K) 

~TN A . 6 TOJ(K)f{u>{K),T) ]T V g (uf) 

ui'<w(K) 

= TN A .ehW{K)f{w{K),T)Ne{K). (14) 

In going to the second expression on the right hand side the a//(a/,T) is for simplification 
replaced by its intermediate value at ZJ « ^(0 + ui), on the grounds that in the narrow range 
(0, w(Ki)), u)f(uj,T) varies insignificantly. The E p h,g(ui) above is equal to the loss of E s 
over the distance £: E s {£) — E s (0). Now, the average viscous force, Fg, is determined by 
Fg = E p h.g/£, where the E p h.g as defined in Eq. ( [LI] ) in turn equals the total dissipative 
work EdQ^EdSl done by the wall to the flow. Next, by its definition, and by combining 
with Eq. (UJ), the viscosity between the wall and the flow is then obtained as 

Fg _ Eph.g 



Thu;(v s )f(u;(v s ))Mg(u;) 



(15) 



where J 21 (~ '" s ~° ) is the lateral velocity gradient and Ag = N A .g ■ a 2 the area of the flow- 
wall interfacial layer. Furthermore, superfluidity experiments have commonly shown that 
an unconfined superfluid exhibits essentially a normal viscosity which may be represented 
here by that of He I, denoted by 7/Hei- Thus, 770 ~ ?7Hci- Comparing the equations ( |l5| ) for 
8 = 2 and respectively for a fixed £ and placing in them the preceding relations, we then 
obtain the relative viscosity of a 2D SC superfluid flow, namely the ratio of 772 to 7y HoI , to be 



7? Hel ~ ryo M a {K) 2ir \L/d) 



(16) 



In going from the second to the third expression on the right hand side above, the parameters 

Sv 
8x ' 



r, £, p-, a 2 , and f(to,T) common to the samples of 9 = 2 and 9 = cancel out, yielding a 
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r2,o as of Eq. (ill). This is to say that, the extraordinary reduction of the number of states 



of a confined superfluid as observed in Sec III A in effect accounts for the reduction in its 
viscosity For L/d « 10 3 - 5 • 10 4 : 

< 10- 7 ~ 1(T 10 , ©' 

which represents for the given dimension of the confined He II sample the lowest relative 
viscosity namely that due to the pure superfluid. In the superfluidity experiments where 
the superfluid sample used has a dimension ratio L/d ~ 5 TO 4 , it has been estimated [§M[l9l 
that the viscosity of the superfluid is lower by at least a factor 10 -6 compared to that of a 
normal fluid, or, that of He I; i.e. % XP / 7 ?Hci — 10~ 6 (circle in Fig. |^a). In comparison, our 
theoretical value (solid or dashed line in Fig. 0a) is well below this. Since a small portion 
of the normal fluid component would contribute to the experimentally measured viscosity, a 
somewhat lower theoretical value for the pure superfluid as obtained above is as expected. 
Therefore, we may conclude that the theoretical 772 reproduces the characteristic order of 
magnitude of experimental data. 

For a ID confined superfluid, the relations corresponding to Eqs. (0) and flll] ) at K = K\ 

are: M\{K\) = ^ + 3, and r\$ — & (x7d) respectively. This says that, similar to Afa(Ki), 

the reduction of Ai(-fCi) is, too, proportional to (d/L) 2 ; and 7*1,0 (dashed line in Fig. ||a) is 
only twice larger than 7*2, Oi which makes no essential difference to the absolute scale of r^o- 
In the case of 3D SC, Nz(K < K\) = and hence 773 = 0. Any measurable viscosity should 
therefore result from the normal fluid component of He II. The 3D SC may be compared to 
the superfluidity experiments where He II flows through porous materials [ §,|o|; although 
experimental viscosity of He II (below v c ) so far given is not well separated as to resulting 
from which component. It is suggestive however to observe that in the experiment by Mehl 
et al.[ [| who employed a 3D SC, the ratio i]T P / V^ci was estimatcd[ |l9| to be at least 10~ n , 
which is an order of 10~ 5 better than the estimated ratio T^/i/Hei < 10~ 6 for the capillary 
experiments [ 0,0] with a similar d value. 

The theoretical evaluation of viscosity above, supported by its good agreement with ex- 
perimental results, has already underlined a superfluidity mechanism of 4 He which we now 
summarize as follows: In a superfluid subject to a 2D or ID spatial confinement with a 
dimension ratio d/L, the number of phonon states ]\[q(K) (9 = 2 or 1) will be reduced to 
the fraction (d/L) 2 of that of a unconfined sample, and is effectively zero for d/L << 1. 
Whereas for a 3D SC, Afs(K) = 0. In any of the cases, the flow motion of the superfluid 
with e s satisfying Eq. (j^j will be effectively "non-dissipative" on the time scale during which 
a normal fluid flow velocity will have long been exhausted. This mechanism has its origin 
in Quantum Confinement Effect (QCE), and will be referred to as the QCE superfluidity 
mechanism. 

That the QCE is responsible for the superfluidity motion of He II is further directly pointed 
to, in the first instance, by the experimentation by Reppy et al.[|2l[] and others[ ^,p2[, which 
showed that after being initially set to rotation in a vessel, a He II current persisted typically 
for hours. This clearly reveals that the "persistency" has to be of a "non-dissipative" origin, 
or, the dissipation is negligible at least within a duration of hours. In the second case, a wide 
range of experiments show that He II actually displays normal friction in a macroscopic bulk 
flow, its non-dissipative superfluidity only presents in a narrow channel - direct evidence of 
the confinement effect - and below a critical velocity, v c (d) which is inversely proportional 
to d (circles in Fig. ||). The good quantitative agreement between the theoretical value of 



v c (d) to be obtained in Sec. IV B and experimental data will add one further support to 
the QCE superfluidity mechanism. Indeed, it is also as expected from basic principles of 
quantum mechanics, that QCE is deemed to play a role in the superfluidity phenomenon, 
which does occur in liquid helium when it turns into a quantum liquid, the He II, and is 
spatially confined. 



IV. CRITICAL VELOCITY Vc- 
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A. Prediction of the presence of a v c 



When the e s of Eq. (|j) begins to just exceed E Q A ph , the states satisfying K > K x will be 
progressively accessible to excitations. Make reference to the geometric relation indicated 



i = (nK z y ~ (iK x y 



and n = with i = 1, j = and v = n — ijj- = 1 in the plot. Similarly, a grid line of (i, j) 



by the large triangle OAB (upper graph of Fig. [Lja), where ABi = \ — 

d 

> (1,1) (not shown in the figure) has the length CD 

straightforward that the A/2 (if) and T>2(u>) within K < K n = iK x + vK z = nK z (where 
v = 1, 2, . . . , ^) can be obtained to be 



nK z ) - (V^+fK x ) . It i 



V — — n 



<y,r,v) 

= (0,0,1) 



[4j + 8 E ABjd] + [E g V j. (4 + SCDjj)] 



K z 



[2n + 1] 

z —i 

£ 

i' =i,i'>i' 

„>V.' 2 +j' 2 i 



4iU = ,, # +8 £ 



and 



4 + 8^/n 2 - ( I 



, 9n . 



P 2 ( w (if))/(— ) = 2 + 4| n=i + £ 



8?; 



8n 



>' =i,f>j', 



+ 4 5 



(&&) 



where ii^i < X < (i + l)Ki, 



lb) 



and <^ = £,f,§,2,2,3, for = (1,1), (2,1), (2,2), (3,1), (3,2), (3,3). For i = 1, the 
above simplify to: 



and 



M(A') = (2n+l)+4-<| i= i+i 



0, 



K X <K< y/2Ki 

4 • z| l=1 + 8,/« 2 - (V2§) 2 , y/2Ki<K<2K u 



0. 



S/i 



V« 2 -(V2^) 2 ' 



V 2 (uj)/(dn/duj) = 6 + 



K X <K< splKx 
V2Ki < K < 2Ki. 



8n 



-(f) 



£Y2 



(8.6)' 
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For l/(d/L) = 5 x 10 4 , the resulting Af2(K) and T>i{ijS) (g„/g^)jv are graphically shown in 
Fig. §>|c, where N a = ^ = 1.07 X 10 15 . N 2 {K) is a stepwise function (Fig. |.b); it 
first increases discontinuously at K = K x by a small amount 5N 2 {K X ) = 4, and then at 
K = K x + 1 • K z by a large amount: 

SM 2 (K 1 + K z )= Af 2 (K x + K z ) - N 2 {K X ) 
L/d+l r— 

- E Sw = aJ T + i. (17) 

l'=0 

i.e. 5J\f 2 (Ki) amplifies in proportion to y/L/d; K z , being << K x , is dropped from 8J\f 2 . For 
the given L/d, 8N 2 (Ki) is as large as ~ 2.5 x 10 3 per atom (inset of Fig. ||b). After this 
stepwise turning, J\f 2 (K > Ki+K z ) rises much more rapidly with K as compared to Af 2 (K < 
K x ). Similarly at each higher K xl = i(±)K z (^'V ') or K n = ^Tf(±)K z = 

(1,1),.. .)value, M 2 (K) has a discrete jump by 8N 2 (iK x ). These above discontinuous jumps 
of M 2 {K = iK x ) correspond to sharp peaks in the T> 2 {uj(K)) v.s. n function (Fig. ||c). 
Similarly, for a ID SC at K = nK z and i = 1: N\(K) = 5o >v + Y^,'\v~oa ^i'i'V = 



can be then similarly obtained. It can be noted that M 2 {K > Kx) 
or Afi(K > Ki) is still extremely small when compared to JVq(K > Ki). 

From Eq. (J15|), it directly follows that r\ 2 will behave analogously as Af 2 (K) described 
above. Such behaviour indeed resembles that of thermal excitation related quantities as 
observed experimentally, including pressure difference [ p3|-p5[, or temperature difference [ 
across a flow that is required to produce a small increase of the flow velocity, or both [ 
28y, which have shown to rise rapidly when v s begins to exceed v c . The close resemblance 



26 



27 



is illustrated by plotting in Fig. |2jb the experimental overshoot data (the amount of He II 
flux driven by the temperature difference across the channel generated by a previous heat 
supply) from Ref ^ (circles) together with the theoretical Af 2 (K) indicated by the solid 
curve. If suppose that now a constant pressure difference is maintained to increase the flow 
velocity. Then, by (@.b), when v s exceeds a v c corresponding to K x the flow will receive 
a sudden deceleration so as to appear as if suddenly disturbed. This actually resembles 



another typical phenomenon observed in experiments [ |l9[. As we indicated in Ref. [ 29 1, a 
direct measurement of the peaked structure in T> 2 (K) (in the low K region) using neutron 
scattering is not presently feasible at least for a macroscopically confined He II, as is limited 
by instrument resolution (5q > 0.005 A -1 and Slo > 0.01 /zeV being roughly the highest 

achievable today) compared to the small K x and A p h involved; and similarly also for a 
microscopically confined He II if lower resolution is used. In the case of the recent neutron 
scattering experiments[ |l5|], where He II is e.g. confined in pores of size ~70 A in Vycor, 
the resultant K x , ^0.09 A -1 , cannot be resolved with the q resolution > 0.4 A -1 used in 

these experiments. Furthermore, although with a high resolution instrument the K x and 
Aph due to pores < 100 A can in principle be resolved, however the possible effect from the 

varying size and shape of the pores in the porous material, which if non-trivial will smear 
the information from each single discrete state, would need be first investigated. 



B. Exact evaluation of the v c (d) 



Consider now the flow velocity is increased such that th e e s of each collision atom is just 
able to produce an excitation of the states at K\\ by Sec. IV A , here v s has its critical value. 
For the critical condition in question here, we only need to consider head-on collisions and, 
of these, the successful excitations. Then, the threshold excitation energy is given by the 
summation of Eq. (Ill]) over only the degenerate modes at K = K\ : 
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(£ph. 2 (ifl)-ffph.2(ifl-^)) 

Sph - 2 TN A . 2 f(ZJ(K),T) 

= hu'V 2 (uj') (18) 

uj' —uj{K\) 

where Y,u'=u{Ki) hu"D g {u') = J2 a A P h = e P h is as in Ec f ©! a = ± ( x ' 2/' z )- 0n the otner 
hand, for a flow atom in collision to be able to produce the minimal excitation above, in the 
severest case the atom has to lose its entire translation energy as given by Eq. (|l2|) and be 
decelerated from v s — v c to v s = 0. Combining Eqs. ([l2]) and (|l8|), we therefore obtain the 
threshold relation for the flow to thermal energy conversion, or the condition for the onset 
of critical velocity v c (= (v s ) m i n ), expressed for per collision atom to be 

0- \m d v 2 c = A ph , or e c s = e c ph . (19) 

a 

Two remarks may be made here. First, the collision atom having a large effective mass 
represents a classical particle; its collision with the wall, being inelastic, is thus governed 
only by the energy conservation condition in Eq. ([l9]). This is in contrast to the un- 
corrected quantum helium atoms in Sec. |y|. Second, the "stopped" atom as in Eq. ( fL9| ) 
will exchange momentum with its surrounding atoms in the flow thereby reducing the total 
flow translation energy; since on satisfying Eq. Jl9| ) many atoms at a time will be "stopped" , 
this will then decelerate the flow motion abruptly. 

Suppose now the channel of a cylindrical shape has a diameter d and length L, and thus 



an interfacial area of A — irdL and volume V = L^-. We thus have Na = = 



mil <uiw Minimi- i — J. W " 1™-™ AT . — A — ndL 

M s = N s m = = ?f£m, and 

M s md 



N A 4a- 




Channel Width d (m) 

FIG. 3. Critical velocity v c {d) versus channel width d for He II. The lines indicate theoretical 
v c (d) values predicted from the QCE superfluidity mechanism for a 2D confined superfluid bulk 
(solid line) and for non-interacting single He atoms of an effective mass (dashed line). Circles 
represent experimental data (at about 1.4 °K) from Ref. [H^. flow, 

Substituting in Eq. (|1]) with the expressions (|l3|.b) and (0)' for m ef and A p h respectively 
and reorganizing, we obtain the critical velocity 



/48h Cl a 1 6.42 x 1CT 7 , , , 

r, =\l — - = m/sec . 20 

m d d 
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The parameterization above, where d in meter, is obtained by using c\ = 239 m/sec, a — 
3.6 x 10~ 10 m (He II interatomic spacing), and m = 6.64 x 10~ 27 kg ( 4 He atomic mass). 
The resultant v c (d) (solid line in Fig. 0), containing no adjustable parameters, agrees 
satisfactorily with the experimental data (circles) for wider channels down to d ~ 10~ 6 
m. We accordingly refer to d < 10 _6 m as narrower or microscopically narrow channels for 
which the v c (d) predicted above deviates from the experimental data, and suggest that the 
deviation has to do with the fact that d becomes comparable to or less than the correlation 
length of the superfluid atoms. We can make a rough check, by using the thermal de Broglic 
wavelength A = /i/V27rmfcsT w 9 A at T — 1 °K as the lower limit of the correlation length, 
due to only atomic wave overlap [ |32]|, ks being the Boltzmann constant. For d = 10~ 8 m, 
for instance, the d encloses in itself as few number of A's as d/A ~ 10. In contrast for a 
wider channel, of e.g. d = 10 -4 m, d/A ~ 1 x 10 6 . 



V. SUPERFLUIDITY IN NARROWER CHANNELS 



In a microscopically narrow channel, the superfluid atomic interaction will reduce with d. 
Accordingly the phonon description of the excitation, which is evidently a result of a strong 
many-atom correlation, will become less adequate; rather, the excitation more appropriately 
involves single atoms. We assume this is the case with d < 1CP 6 m. Within the small d, on 

the other hand, each helium atom will directly sense the channel boundaries and collisions 
will primarily occur between a helium atom and the channel wall. To include the two-fold 
effect above, the atoms may first be approximated as un-correlated but each having an 
effective (heavier) mass, m* , to account for the small but finite correlation. Second, each 
single atom is now confined to d in the lateral directions in real space, and confined to a 
potential well of an infinite depth and a width d in the atomic wavevector k space. The 
single atom standing wave function for this is Vk( r ) = Ck& + C-k&~ ■ The boundary 
conditions are k X i = ik x , and k z i = lk z , i, I = 0, 1, . . .; where k x = § and k z — j-. Applying 

these to the wave function and solving gives the eigen kinetic energy e n — -g^^p, where 
n = 1, 2, . . ., for the lateral k x and k y states. The minimum excitation energy is 

3/? 2 

A a = 6 2 - ei = -^. (21) 
8m* gt 



The threshold excitation condition as of Eq. (|19|) or (|6j) now applies equally to the atoms, 
but with the A p h replaced by the A a . Accordingly, the QCE superfluidity mechanism 

for the single atoms, of a second type here, can be similarly stated as that for the phonons, 
in terms of the reduction of atomic wavevector states. Similar to the preceding discussions, 
a critical velocity can be defined when the flow energy E s = ^N s mv 2 = N s e s begins to 

exceed the (total) threshold excitation energy N s where a = k x , k y , k z so that e c a = 

V,A„ 3A,. T':r.i> 

- mv 2 r = -. 22 

2 c 8m*d 2 v ' 

The critical velocity follows to be 

3h 5.1 x 1(T 8 

v c (d) = —— ~ , (23) 

Zyfjmd d 

where 7 = m*/m. For the parameterization above, we have made the crude estimation: 
7 ~ (A/a) 3 = (/i/aV27rmfc B T) 3 ~ 8.6, with T = 1.4 °K. Values for other parameters are 
as specified earlier. The resulting v c (d) (dashed line in Fig. |3|) is seen to have an improved 
agreement with the experimental data in the low d end. There is however still a noticeable 
discrepancy; it is evident that first of all, m* needs to be more realistically determined. 
Furthermore, as d is reduced to the atomic scale, we expect that the localization of the 
center of mass of each atom wave in d, i.e. that the atomic translation is suppressed to 
vibrational motion, can result in an even smaller v c . 
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VI. AUXILIARY COMMENTS 



Supplementary to the preceding discussions, the following points are noted. 1). The scarcity 
of the K modes obtained in this work has no connection to the "Bose-Einstein condensa- 
tion" . As we have noted, the modes dissatisfying Eq. are both thermally absent and 
inhibited from excitation on perturbation; whilst all of those satisfying Eq. (H) are never- 
theless thermally regularly excited. 2). The well defined discretized excitation state of the 
superfluid, being an experimental fact and having been shown here to lead to superfluid- 
ity, implies a fully coherent atomic motion, which in turn implies relative localization of the 
atoms in a bulk superfluid. In contrast, the state discretization, which may present at a long 
wavelength in any fluid, cannot sustain its effect in an essentially diffusive atomic system, 
as is in the classical fluid or in a bulk quantum fluid which is assumed to be consisting of 
free atoms. 3). We finally comment that the superfluidity criteria as given by Landau and 
by Feynman have both incorrectly neglected the primary low energy excitations of phonons. 
Landau's this ignorance can be traced to the unphysical comparison [^| between velocities of 
different masses - the v s of the correlated superfluid atoms having a mass (N s /Nji)m where 
N s /Na > 10 3 , and the ci of a phonon having a zero mass[ |33f ; in the case where a fluid bulk 
collides with the wall inelastically, such a comparison does not conserve energy . Consistent 
with his conceptual error, Landau's critical velocity (~ 100 m/sec) is at least 10 3 times larger 
than the experimental value and contains no mechanism for d dependence. Feynman in his 
treatment [ (1955) used Landau's critical velocity (vl < Ca, C2 being the phase velocity 
at K w 1.93A -1 ) as a "valid" starting point, and consequently ignored phonon excitation. 
Furthermore, Feynman's attribution of turbulence as the source for breaking-up superflu- 
idity suffers the pathology that, for a superfluid assumed to be " non- viscous" , turbulence 
cannot be produced until the turbulence has already been produced. 
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